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Star junctions and watermelons of pure or random quantum Ising chains : 
finite-size properties of the energy gap at criticality 
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We consider M > 2 pure or random quantum Ising chains of N spins when they are coupled via a 
single star junction at their origins or when they are coupled via two star junctions at the their two 
ends leading to the watermelon geometry. The energy gap is studied via a sequential self-dual real- 
space renormalization procedure that can be explicitly solved in terms of Kesten variables containing 
the initial couplings and and the initial transverse fields. In the pure case at criticality, the gap 
is found to decay as a power-law Am oc with the dynamical exponent z(M) — ^ for the 

single star junction (the case M = 2 corresponds to z = 1 for a single chain with free boundary 
conditions) and z{M) = M — 1 for the watermelon (the case M = 2 corresponds to 2 = 1 for a 
single chain with periodic boundary conditions). In the random case at criticality, the gap follows 
the Infinite Disorder Fixed Point scaling In Am = —N'^g with the same activated exponent ip = ^ 
as the single chain corresponding to M = 2, and where g is an 0(1) random positive variable, whose 
distribution depends upon the number M of chains and upon the geometry (star or watermelon). 


I. INTRODUCTION 


The quantum Ising chain is the basic model for quantum phase transitions at zero temperature [l| : the pure chain 
is exactly soluble 0 and is the same universality class as the two-dimensional classical Ising model via the classical- 
quantum correspondence [ll ; the random chain is also exactly soluble via the Ma-Dasgupta-Fisher Strong Disorder 
renormalization approach [3| and is the prototype of Infinite Disorder fixed points with unconventional scaling laws 
with respect to the pure case (see 0 for a review). 

Remarkably, for an open Ising chain containing N quantum spins with arbitrary transverse fields h{i) and arbitrary 
couplings J{i) 
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there exists an^ej^licit exact formula for the surface magnetization at one end i = N when the spin at the other end 
z = 1 is fixed 
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For the disordered chain where the transverse fields h{i) and the couplings J(i) are random, the surface magnetization 
thus involves the specific structure of Kesten random variables satisfying simple multiplicative recurrence [iS3- 
Kesten random variables also appear either directly in many other discrete disordered systems [lll - [l^ . or via their 
continuous analogs given by exponential functionals of Brownian motion [20l - l^ . The probability distribution of the 
surface magnetization over disordered samples has thus been analyzed in detail • 

Unfortunately, other interesting observables are not given by explicit exact formulas in terms of the initial disorder 
variables, but various proposals have been made to relate them to the known surface magnetizations : 

(i) for the chain of Eq. [T]with free boundary conditions, the end-to-end correlation 


C™’'/(l,Af) =< CT^(l)a^(iV) > 


( 3 ) 


whose statistics has been computed via the appropriate extension of the Strong Disorder Renormalization Q, 
was found numerically to be very close to the product of the two surface magnetizations Q 


C>™r/(i,jv) ~ m™’’'^(Af|l)m™’’'^(l|iV) 


( 4 ) 


(ii) the energy gap between the ground state and the first excited state of the chain of Eq. [T] with free boundary 
conditions 


An = El — Eq 


( 5 ) 
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whose statistics has also been studied via Strong Disorder Renormalization (23|, was found to be related to the 
end-to-end correlation of Eq. [3] via 


N 


N-1 


In Ajv - In (1, A) ~ ^ In h{k) - ^ In J(fc) 


( 6 ) 






as a consequence of the structure of the Strong Disorder renormalization rules 23|. On the other hand, using the 
free-fermion approach, the gap was found to be related to the product of the two surface magnetizations 
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provided An decays faster than 1/iV [2^. Eq|6]and Eq. [7] are equivalent if one takes into account Eq. |4l 

In this paper, we obtain via real-space renormalization some generalizations of these formula for star junctions of 
M quantum Ising chains, and we use them to analyze the properties of the energy gap. More precisely, we consider 
M quantum Ising chains containing N spins crp{i) each with 1 < i < N and I < p < M 
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when they are coupled via a Star junction involving an additional central spin u'(O) 
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or when they are coupled via two Star junctions involving two additional spins cr(0) and a{N -f- 1) leading to the 
Watermelon geometry 
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The star junction of Eq. |3[with random transverse fields and random couplings has been studied recently via Strong 
Disorder Renormalization |25l | to obtain the disorder-averaged magnetization at the center of the junction cr^(O) in 
the thermodynamic limit of semi-infinite chains N —>■ -|-oo. Pure Quantum Ising chains with uniform transverse field 
hp{i) = h and uniform couplings Jp(i) = J when they are coupled via a Delta-junction 
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have also been studied in the limits of weak coupling Jp_p/ —>• 0 [26|, l27| and of strong coupling Jpy —>■ -Pc xi l26ll. The 
Delta junction of Eq. [TT]is expected to have the same physical properties as the Star junction of Eq. [31 In , the 

case of two Delta junctions at the two ends has been also considered and is expected to be similar to the watermelon 
geometry of Eq. [TU] involving instead two star-junctions. Let us finish this introduction by mentioning that other 
quantum models have been studied on the star-junction geometry [2^, [2^. 

The paper is organized as follows. In section [iTl we introduce a real-space renormalization procedure for the 
star junction that can be explicitly solved in terms of the surface magnetizations of Eq. [3] and other related Kesten 
variables, and we obtain generalizations of Eq. |3]for the correlation between the end spins, and of Eq. |7]for the energy 
gap. In section imi we describe the results for the pure star junction, and in particular the dynamical exponent z{M) 
at criticality. In section IIVI we discuss the results for the disordered star junction governed by an Infinite Disorder 
Fixed point. Section |V] is devoted to the watermelon geometry. Our conclusions are summarized in section lyTl 


II. SEQUENTIAL SELF-DUAL REAL-SPACE RENORMALIZATION 

The Fernandez-Pacheco self-dual real-space renormalization introduced initially for the pure quantum Ising chain 
[s^ has been recently extended to the disordered chain (sil - l^ . As in any block renormalization procedure, the 
renormalization rules are applied in parallel to all blocks of size 6 = 2. Here we propose to apply the Fernandez- 
Pacheco self-dual renormalization rules sequentially around the central site of the junction of Eq. [31 
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A. Renormalized Hamiltonian 
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Figure 1: Illustration of the renormalization procedure for the star-junction of M = 3 chains of length N : 

Top : Initial Hamiltonian of Eq. [9] with the transverse fields hp{i) and the couplings Jp{i) 

Bottom : Renormalized Hamiltonian of Eq. [T^] after n elementary RG steps, where the renormalized center is characterized by 
the renormalized transverse field (0) and the renormalized couplings (0). 


More precisely, after n RG steps, we consider that the initial Hamiltonian of Eq. IH] has been renormalized into (see 
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Fig. [T|) 
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where all the spins n+l<i<A'^of the M chains have kept their initial transverse fields hp{i) and initial couplings 
Jp{i) of Eq. |51 and where the central renormalized spin a{R^{0)) with its renormalized transverse field (0) and 
its renormalized couplings (0) replaces the central (1 + Mn) initial spins, namely the initial central spin and the 
initial first n spins of the M chains. The initial Hamiltonian of Eq. IH] corresponds to n = 0 so that the initial condition 
of the RG flow reads 

/i«”(0) =/i(0) 

jf(0) =Jp(0) (13) 


B. Renormalization rules 


To go from to , we decompose into the sum of two contributions 
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The Hamiltonian is diagonal in CT^(i?”“^(0)), and for each value g^(R"~^(0)) = ±1, one has to diagonalize 

independently the single-spin effective Hamiltonian for CTp(n). We refer to for more details and simply state the 
final output : the projection of the complementary Hamiltonian V onto the two degenerate ground states of 
leads to the following RG rules [s^ for the central transverse field 


and for its couplings 



(15) 


(16) 


C. Solution for the renormalized conplings 

It is convenient to rewrite Eg 1161 as 
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in order to obtain the simple iteration 
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so that the final result reads 


Jp«"(0) 


Jp{p^^ 
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^p(fc) 
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The parenthesis on the right hand side coincides with the surface magnetization m®“’’^(n) at z = n of the p chain 
when the magnetization at the center z = 0 is fixed (Eq. [2]) 


m;“’'^(n|0) 

So the solution obtained in Eq. [12] reads 
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( 21 ) 


Note that this formula also suggests that the sequential RG that we have described is actually somewhat equivalent 
after n steps to the large-block self-dual procedure involving a block a size n studied in Ref [34l |. 


D. Solution for the renormalized transverse field 


It is convenient to consider the ratio obtained from Eqs [12] and [T51 
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in order to obtain the simple iteration 
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The solution obtained in Eq |2T] for the renormalized couplings in terms of the surface magnetizations of Eq. 
yields for the renormalized transverse field 
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Using the explicit form of the surface magnetization of Eq. [20l this can be also rewritten as 


M 


‘(o)=Mo)n 


hp(l)hp{2)...hp(n) surf 


';J-^\jpiO)Jp{l)...Jp{n-l) 


(n|0) 



/zp (l)/ip(2). ..hp (n) 
Jp(0)Jp(l)...Jp(n-l) 


1 + 


hl{n) hlin - 1 ) 

Jp{n — 1) Jpin — 1) Jp{n — 2) 


hl{n) hl{l) 


J2(n-1)-J2(0) 


^p(l)^p(2)-/ipH hl{l)hl{2)...hl{n - 1) 


4(0) , . 


p=l . 


1 + y 

hL\ 4(k) 


(25) 


E. Gap of the star junction with free boundary conditions 


For the finite star junction of Eq. [9l the renormalization procedure described above ends after N RG steps where 
the renormalized Hamiltonian of Eq. 1121 contains a single renormalized spin a{R^{0)) representing the whole star 
junction 


Hp^ =-h^ {0)a^{R^{0)) 


(26) 


The remaining renormalized transverse field (0) corresponds to the energy scale associated to the global flip of 
the junction and thus represents the energy gap (up to an unimportant factor of 2) 


Ap^ = h^ (0) 


Eq. [24] yields 
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(28) 


This formula is thus the analog of Eq. |7| : the numerator contains the product of the M surface magnetizations and 
the product of all the transverse fields, whereas the denominator involves the product of all the couplings. 

The alternative form of Eq. [25] reads 
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F. Correlation between the end spins 

Before the last Hamiltonian of Eq. the renormalized Hamiltonian at step n = N — \ (Eq. IT^ 

M M 

Rp-r_, = -h^^-\^)a\R^-\p) -Y^jf'\^)a\R^-\p)al{N) -Y,hp{N)a;{N) (30) 

p—1 p—1 

involves the initial end-spins <Jp{N) of the M chains, and the renormalized central spin (T(i?^“^(0)) representing the 
rest of the sample. The eigenfunction corresponding to these two last RG steps reads (see [Hj for more details on the 
eigenfunctions of of Eq. [HI) 
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The joint distribution of the M end-spins Sp{N) is thus obtained after the integration of the renormalized central 
spin 


V{Sp{N))= ^ ,p\So-,Sp{N)) 
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(32) 
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In particular, the correlation between the M end-spins reduces to 
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and can be thus rewritten using Eqs [201 and HI] as the product of the M surface magnetizations when the center is 
fixed 


M M 

<n5p(A^)>=nK“'’^(^io) (34) 

p—i p—i 

This formula is thus somewhat the analog of Eq. [H 

III. RESULTS FOR THE PURE JUNCTION OF M > 2 CHAINS 

For the pure case, the control parameter of a single chain is the ratio between the ferromagnetic coupling J and 
the transverse field h 



The exact surface magnetization of Eq. [2] reads 


(35) 


^W(jv|l) 
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A. Ferromagnetic phase A > 1 

In the ferromagnetic phase A > 1, the surface-magnetization of Eq. [36] converges towards the finite limit 

^ -Poo) = [1 - A-2] 5 (Xk^i+ {K - I)'^'’ with I3s = ll2 
As a consequence the renormalized coupling of Eq. [21] also converges towards a finite limit 

jE( 0) = Jm"“’'^(A|0) oc 

N^+ao ■' 

whereas the renormalized transverse field of Eq. HS] converges exponentially towards zero as 

(0) = h 


A2- 1 
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B. Paramagnetic phase K < 1 


In the paramagnetic phase K < 1, the surface magnetization of Eq. 1361 decays exponentially with the distance 

K-^ - 1 




K--^n _ I 


OC = 

N—¥-\-00 
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and involves the correlation length exponent v = 1. As a consequence, the renormalized coupling of Eq. 1211 behaves 
similarly 


J«"(0) = OC 

^ N—>-+oo 

whereas the renormalized transverse held of Eq. converges to the hnite limit 

M 
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C. Critical point = 1 

At the critical point = 1, the surface magnetization of Eq. [551 decays as the power-law 

= N-i (43) 

As a consequence, the renormalized coupling of Eq. decays with the same exponent 

J^'^(O) oiN-^ (44) 

whereas the renormalized transverse held of Eq. representing the gap decays as the power-law 

Aw = hf^=i(0)ociV-^W (45) 

with the dynamical exponent 

z{M>2) = ^ (46) 

For M = 2 corresponding to a single open chain with free boundary conditions, this result reduces to z{M = 2) = 1 
as it should. For M > 2, Eq. [46] is in agreement with the strong-coupling perturbative expansion of Eq. (12) of 
Reference 


IV. RESULTS FOR THE DISORDERED JUNCTION OF M > 2 CHAINS 

For the disordered junction, the usual control parameter Q 

^ _ In h — In J 

var(hi h) + variln J) 

valid for Gaussian distributions can be more precisely dehned as the root of the equation [ssl - ls^ 


(47) 




(48) 


that makes the link with the literature on Kesten variables imi. It is also usual in the held to redehne the length 
scale to absorb the variance Q, in order to have the single parameter 5. 
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A. Reminder on the statistical properties of the surface magnetization 

Let us briefly summarize the known results on the statistics of the surface magnetization of random chains [3, 
in the universal critical region around d = 0 : 

(1) In the paramagnetic phase <5 > 0, the surface magnetization reads [ 3 , 13 | 

In (N) = —2SN + yfNv 

where u is a zero-mean random variable distributed with the Gaussian distribution 

1 2 

G{v) = 

In particular, the typical value 




e 


lnm*“’'/(iV) = -2bN 


(49) 

(50) 

(51) 


involves the typical correlation length exponent vty^ = 1, whereas the finite-size scaling properties determined by the 
balance between the two terms of Eq. |49]are governed by ^s_= 2. 

(li) at criticality (5 = 0, the surface magnetization reads [3, 

(N) = —'/Nw 

where w is a positive random variable distributed with the half-Gaussian distribution 

e{w > 0 ) 

F[w) = - -j =—e 

VTT 

(iii) in the ferromagnetic phase d < 0, the surface magnetization reads [3,0 

Inm^"’'/ {N) = -B 

where is a positive random variable distributed with the exponential distribution 

Q{B) = 9{B > with ^ = —2S 


(52) 


(53) 


(54) 


(55) 


From these known results concerning the surface magnetization of a single random chain, we may now derive the 
statistical properties of the renormalized parameters for the random star junction. 


B. Random star junction at the critical point <5 = 0 


At the critical point i5 = 0, the renormalized coupling behaves as the surface magnetization of Eq. 

in J^Q oc In (N) = — Vnw 

The gap Am{N) given by the renormalized transverse field will be also of the form 

In Am(A^) = InhfJo = -y/Ng 
and thus involves the same activated exponent 

^ = l 


(56) 

(57) 

(58) 


as the chain corresponding corresponding to the special case M = 2. The variable g = X)p=i Eq. [57] is 

distributed with the convolution of M half-Gaussian distributions (Eq. (531) 
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(59) 
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However this expression is expected to be valid only for sufficiently lar ge g , whereas the region of small g requires a 
more precise analysis, as shown by the Strong Disorder calculations of concerning the case M = 2 of the chain. 
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C. Random star junction in the Ferromagnetic phase <5 < 0 


In the ferromagnetic phase (5 < 0, the renormalized coupling remains finite as the surface magnetization of 
Eq. [54l 


In oc lnm^“’'/ (TV) = -B 
whereas the renormalized transverse field will behave as 

In = —2SMN + 'J MNv 

where u is a Gaussian random variable distributed with Eq. 1501 


(60) 

(61) 


D. Random star junction in the Paramagnetic phase 5 > 0 

In the paramagnetic phase d > 0, the renormalized coupling inherits from the behavior of Eq. 
magnetization 

In oc In (N) = —26N + y/Nv 
whereas the renormalized transverse field remains finite 

ln/if>o = -B 

where B = Bp is distributed with the convolution of M exponential distributions of Eq. [55] 
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I of the surface 
(62) 

(63) 

(64) 


V. STUDY OF THE WATERMELON GEOMETRY 

In this section, we use the results of the previous sections concerning the single star junction with free boundary 
conditions (Eq. jH]) to analyze the watermelon geometry of Eq. |Tni where the M chains are coupled by two star 
junctions at the two ends. 

A. Real-space renormalization for the Watermelon geometry 

We apply the same renormalization rules as in section [nj around each of the two junctions, so that after ^ RG 
steps, the renormalized Hamiltonian reads 


Tjwater melon 
£1 N 

R~T 


= (0)a"(i?-(0)) - (TV + l)a^{R-{N + 1)) - a%R-i0))a^ {R-{N + 1)) (65) 


The renormalized transverse fields of the two centers of the junctions are obtained from Eq. 
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J,^(0)J,^(1)...J,^(^-1) J,^(0)J,^(l)...J,^(f-2) 
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11 


and 
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whereas their coupling is obtained by summing the contribution of the M chains 
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surf 
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Jlik) 


1 + 
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'N 


J2(N+1)J2(M+2)-JI{N) 


(70) 


B. Solution for the pure watermelon 


For the pure case, where the surface magnetization is given by Eq. [36] in terms of the control parameter K = J/h 
of Eq. |3S1 the renormalized coupling of Eq. |B5| reads 






w 


= MJ- 


l-K- 


1_^_2-7V 

whereas the transverse fields of Eqs [521 and [571 have the same value 


(71) 
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In the ferromagnetic phase > 1, the renormalized coupling of Eq. 1711 converges towards the finite limit 
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N^+oo 


and the renormalized transverse field of Eq. ^2\ converges exponentially towards zero as 


(0) ~ hn-K-^)^K 
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In the paramagnetic phase K < the renormalized coupling of Eq. [7T] decays exponentially with the distance 

K 

~ MJ{1-K‘^)K^ (75) 

N—^-\-oo 

and the renormalized transverse field of Eq. [25] converges to the finite limit 

N JV M 

(0) = {N+1) =h{l-K^)^ (76) 

At the critical point Kc = 1, one obtains the following power-laws for the renormalized coupling of Eq. [71] 

= MJ-^ (77) 

and for the renormalized transverse field of Eq. 1251 

M 

N N / 1 \ ^ 

h^^{0) = h^^{N + l) (78) 


For M = 2 corresponding to a single chain with periodic boundary conditions, both involve the same exponent, 
that represents the dynamical exponent z{M = 2) = 1 as it should. However for M > 2, the coupling of Eq. |77| 
dominates over the transverse fields of Eq. [751 so that the last decimation corresponds to the projection onto the two 

ferromagnetic states (++) and (-) defining the final renormalized spin representing the whole watermelon, and its 

renormalized transverse field 


represents the gap 


'last 
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’JW 


h 

M 
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(79) 


(80) 


with the dynamical exponent 


zw{M > 2) = M - 1 


(81) 


For M > 2, the dynamical exponent of Eq. [51] is bigger than the dynamical exponent of a single star junction of Eq. 

SHI 


C. Solution for the disordered watermelon 

The statistical properties of the renormalized transverse fields of Eos 1661 and 1671 are the same as in the star junction 
case given in Eqs |S71 [HT] and [HHI with the replacement ^ • 

The statistical properties of the renormalized coupling of Eq. jSSj can be obtained from the properties of the surface 
magnetization recalled in section IIV Al In the paramagnetic phase <5 > 0, it decays as 

N_ _ 

InJ^" = -2SN + VNv (82) 

where v is the minimum value among M variables drawn with the Gaussian distribution of Eq. 1501 In the ferromagnetic 
phase d < 0, it remains finite 

\nJ^^=-B (83) 

where B is minimum among M positives variables drawn with the convolution Q * Q = yF'Be~^^ of the exponential 
distribution of Eq. jSS] At criticality i5 = 0, it decays as 

N _ 

InJ^^ = —'/Nw (84) 

where w is the minimum value among M positives variables drawn with the half-Gaussian distribution of Eq. 1531 
In particular, the gap of the watermelon involves the same activated exponent ^/; = 1/2 and it is only the scaling 
distribution of the gap that depends upon M. 
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VI. CONCLUSION 

For M > 2 pure or random quantum Ising chains of N spins coupled via a single star junction or via two star 
junctions at the two ends, we have introduced a sequential self-dual real-space renormalization procedure that can be 
explicitly solved in terms of Kesten variables containing the initial couplings and and the initial transverse fields. We 
have described the properties of the renormalized couplings and renormalized transverse fields in the ferromagnetic 
phase, in the paramagnetic phase and at the critical point. In particular, we have discussed the finite-size properties 
of the energy gap Am at criticality : 

(i) In the pure case, the dynamical exponent z governing the power-law decay of the gap Am oc was 

found to be z{M) = for the single star junction and z{M) = M — 1 for the watermelon geometry (for the case 
M = 2 corresponding to a single chain with open or periodic boundary conditions, both dynamical exponents reduce 
to z{M = 2) = I as it should). 

(ii) In the random case, the gap was found to follow the Infinite Disorder Fixed Point scaling In Am = —N'^g with 
the same activated exponent ip = ^ as the single chain corresponding to M = 2, so that the number M of chains and 
the geometry (star or watermelon) only enter in the probability distribution of the 0(1) random positive variable g. 

As a final remark, let us mention some link with the recent studies concerning the Shannon and Renyi entropies 
of the ground state wavefunction of the pure quantum Ising chain [H, . The Renyi entropy of index q of the 

quantum chain amounts to study the binding of the ’Ising book’ (see Fig. 4 of 0) where M = 2q half-planes of 
classical two-dimensional Ising models are glued together. Taking the anisotropic limit in the other direction leads to 
the pure quantum star junction discussed in section IIIII Although the present study does not deal with observables 
that are directly relevant for the Renyi entropy, we feel that the expansion of the dynamical exponent z{M) = ^ = q 
for q = l + e around the Ising plane q = t 

z{M = 2{l + e)) = l + e (85) 

suggests that the appropriate finite-size scaling variable could be L'’ = e“ with 

u = eh\ L = [q — l)h\ L (86) 

For the periodic quantum chain, Gregoire Misguich and Jean-Marie Stephan have replotted their numerical data 
showned on Fig. 5 of their paper [s^ in terms of the variable u = (g — l)lnL and have concluded that this finite 
size scaling could be also compatible with their data [l^. For the open quantum chain, the double logarithm (InL)^ 
found numerically at criticality for g = 1 by Jean-Marie Stephan and the singular logarithmic terms in {cst)q^^ 
found at criticality for g ^ 1 via Conformal Field Theory in Ref [d^ also points towards the finite size scaling variable 
of Eq. ES 
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